#ifndef RTNORMCHPP
#define RTNORMCHPP
//  Example for using rtnorm
//  
//  Copyright (C) 2012 Guillaume Dollé, Vincent Mazet (LSIIT, CNRS/Université de Strasbourg)
//  Licence: GNU General Public License Version 2
//  see http://www.gnu.org/licenses/old-licenses/gpl-2.0.txt
//
//  Depends: LibGSL
//  OS: Unix based system


#include <iostream>

#include "rtnorm.hpp"
#include "../../Random.hpp"


class rtNC
{

	public:
		rtNC()
		{	
			int seed = RandomG::getSeed();    
			//T = gsl_rng_ranlxd1;  // generator with good math properties
			//T = gsl_rng_taus; //default generator
			T = gsl_rng_mrg;
			r = gsl_rng_alloc(T);    
			gsl_rng_set (r, seed); 	  

		}

		double operator()(double a, double b, double mu, double s)
		{
    			return VM::rtnormVM(r,a,b,mu,s);

		}
	private:
	const gsl_rng_type *T;
	gsl_rng *r;
  
};

#endif
